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Quantum site percolation as a limiting case of binary alloy is studied numerically in 2D within the tight- 
binding model. We address the transport properties in all regimes - ballistic, diffusive (metallic), localized 
and crossover between the latter two. Special attention is given to the region close to the conduction band 
center, but even there the Anderson localization persists, without signs of metal - insulator transition. We 
found standard localization for sufficiently large samples. For smaller systems, novel partial quantization of 
Landauer conductances, i. e. most values close to small integers in arbitrary units is observed at band center. 
The crossover types of conductance distributions (outside the band center) are found to be similar to systems 
with corrugated surfaces. Universal conductance fluctuations in metallic regime are shown to approach the 
known, theoretically predicted value. The resonances in localized regime are Pendry necklaces. We tested 
Pendry's conjecture on the probability of such rare conducting samples and it proved consistent with our 
numerical results. 
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1 Introduction 

The quantum percolation model exhibits Anderson localization and it was reviewed e. g. in Ref. [1J. We 
believe that it belongs to the same universality class as the Anderson model, thus it undergoes no metal 
- insulator transition (MIT) in 2D, i. e. all states become localized for any positive disorder [2 3] and 
large enough samples. We will regard square, or generally hypercubic disordered samples with L d sites, 
coupled to perfect semi-infinite leads, d is dimensionality. Two characteristic lengths can be introduced, 
the mean free path / and the localization length £. For finite samples and small disorder, with system 
size L significantly smaller than both I and £ we get the ballistic regime and roughly for I < L < £ the 
(quasi)metallic one, with mean Landauer conductance (g) oc and (g) oc L d ~ 2 , respectively. In the 

localized regime L > £ and (g) oc exp(— L /£) . At MIT point, characterized by critical amount of disorder, 
the localization length £ diverges and it remains infinite in full metal. Such MIT surely exists for d > 2. 
In quasi-metallic regime £ is large, but finite, typically for d < 2. The absence of MIT is known also in 
2D Anderson model with orthogonal symmetry, i. e. neither with magnetic field nor spin-orbit interaction. 
The same was numerically checked again [4| for 2D quantum percolation. On contrary, in Ref. [5| and 
quite recently in a paper [6 |, MIT presence at 2D quantum percolation was reported. We aim to support 
the former statement that no indications of MIT are found when calculating the dependence of Landauer 
conductance on system size and other model parameters numerically. 

Concerning experiments, in a recent paper [7 1 the authors report direct confirmation of MIT for 3D An- 
derson model, including correct critical exponent, characterizing the divergence of £. These measurements 
were performed on atomic matter waves, cooled Cs atoms kicked by laser pulses. Such experiment rules 
out problems of decoherence or many-body interactions, thus it is well described by our model. 
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However with real mesoscopic samples, experimenters often observe the MIT on 2D structures, see e. g. 
the review [ 8 ] . But even in our non-interacting, or one-electron models, switching the magnetic field or spin 
- orbit interaction already induces the MIT in 2D. At least the metallic regime seemed to be well understood 
and widely accepted. However, later measurements [9| on quasi- ID gold wires showed, than even in this 
regime the third cumulant of g was non-zero, thus the distribution of conductances differs from the expected 
Gaussian. This can be explained e. g. by electron - electron interaction. One-particle models predict single 
parameter scaling and Gaussian distributions within the conduction band, though with possible exception 
of some isolated points (energies). Alternatively we can have to do with decoherence by inelastic scattering. 
Metallic regime was also measured on 2D electron systems in Si MOSFETs [10], two parameter scaling 
involved, and it was again explained by above mentioned e-e interactions. Thus the single parameter 
scaling theory predictions, especially no MIT in orthogonal 2D systems, are disputed from various sides. 
The non-satisfactory situation was summarized in (8), stating that there is no simple model capable to 
explain all experimental data at once. 

Superlocalization, i. e. over-exponential decay of g(L) oc exp — k > 1, was shown in Ref. ifTTI 
on the percolation cluster at p = p c and explained as a consequence of the cluster's fractal nature. Here p 
is the ratio of randomly removed sites and p c is the percolation threshold. We do not reach the classical 2D 
site percolation threshold p c — 0.407. 

The charge careers, represented by planar quantum waves, have the energy E. The conduction band 
for simple hypercubic sample with the volume L d is the interval —2d<E<2d. Anomalous behavior, 
meaning non-exponential decay of g(L), was reported [12] mainly at band center, i. e. for small energy \E\. 
In Ref. lfl2l Thouless conductivity a was calculated from the dependence of energy spectra on boundary 
conditions and a transition between two types of localization was reported even in 2D. One expected type 
should be sublocalization or weak localization, i. e. power-law decay of a(L) with negative exponent. We 
will test such possible behavior for Landauer conductance g and the largest L available. 

In localized regime, the energy dependences of conductance have resonances, i. e. sharp peaks. They 
ofted reach to conducting values, i. e. g max « 2. The nature of wave functions at resonant energies lfl3l 
has shown that they are of Pendry necklace type, i. e. a chain of comparable localized wave packets. 
Such a quasi ID chain enables the transmission across the sample also in higher dimensional case. Pendry 
predicted [14], that even for d > 1 the probability r of these rare states should decay with L as r oc 
exp(— -\/L/£), a relation he derived analytically in ID. 




As the main tasks, in this paper we aim to look for signs of MIT or sublocalization in the whole con- 
duction band, to test Pendry's conjecture on rare events and to present the properties of the quasi-metallic 
regime, compared even with quantitative predictions. 

2 Model and Method 

Let us now introduce the tight-binding Hamiltonian: 



where i numbers the L x L lattice sites and the second sum goes over nearest neighbors. Later we set the 
hopping V = 1 in arbitrary units, thus fixing the energy scale. With box distribution of random on-site 
energies P(ei), —W/2 < < W/2, the Hamiltonian (Q]) becomes the Anderson model and W is called 
disorder. For a binary system the distribution of the energies has the form: 



where p is the ratio of insulating sites B and it is the measure of disorder in this model. S(x) is Dirac's 
delta function. We can analyze any combination of energies tA and es by our method, but in this paper 





(1) 



<i,j> 



P(ei) = (1 - p)S(e t - e A ) + p5(e % - e B ), 



(2) 
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we concentrate on the quantum percolation model. Hence we set the on-site energy for A-type atoms to 
€a — and the one for the randomly positioned B-type atoms to a large number, mostly es = 10 7 . For a 
completely removed, inaccessible site this energy would be infinite. The difference will be discussed in the 
next chapter. We chose this approach in order to apply the transfer matrix method to calculate the transport 
properties. When we study the eigenenergies, namely their density of states (DOS) p, we simply set all 
on-site energies to zero and also the hopping elements Vy — for all nearest neighbors j of the excluded 
sites i. 

We use standard transfer matrix method for numerical calculations. The square sample is connected to 
two perfect semi-infinite leads. A set of L (generally planar waves enters the square sample from 

the left. They have different wavenumbers in transfer direction k n , n = 1, L, given implicitly by the 
dispersion relation 



E = ~2Vcosk n a-2Vcos- -, < k n < n/a. (3) 

Li -\- 1 

We set the lattice constant a = 1, thus fixing the length scale. The last term is given by hard wall 
boundary conditions in perpendicular direction. The real solutions for k n are called open channels and 
they are possible only for \E\ < 4; the complex solutions lead to exponentially damped waves - closed 
channels and they are omitted at the end. The part of the waves arriving on the right hand side of the 
sample is described by the transmission matrix t, the rest was reflected to the left and it is described by 
the reflection matrix r, both L x L. The elements of the transmission matrix t a p are the amplitudes of 
charge careers transmitting from channel a (on the left) to the channel f3 on the right. The conservation 
law has the form tH + r^r =1, where 1 is identity matrix representing incoming waves and this is a 
generalized Kirchhoff law. The conductance g in units e 2 /h is calculated making use of Landauer formula 
g = 2 Tr tH, the trace restricted to open channels. The 2L x 2L transfer matrix T relates waves on the 
right: vector V r of L outgoing transmitted waves and L (zeros) incoming ones, to those on the left: Vi of 
L incoming and L reflected waves: V r — TVi. It can be shown to fulfill lfl5ll 



T " (r't" 1 tt - rtt-if) (4) 

The main advantage of this formulation is the possibility to cut the sample into slabs, perpendicular 
to the transfer direction. The L transfer matrices of these slabs are then multiplied to get the complete 
T, because outgoing waves of one slab are the incoming ones of the next slab. But we have to do with 
ill-defined matrices, i. e. their eigenvalues differ by many orders of magnitude for samples even much 
smaller than it was the case for Anderson model. In other words, numerical errors, given mainly by closed 
channels, rise quickly with sample size. An efficient way to treat this problem was described in Ref. 1T51 
and at the end of Ref. fl6l . It is based on the fact, that the upper matrices involved in Eq. (O have large 
eigenvalues dominating t , but we need accurate lower eigenvalues to recover t. Thus we multiply the 
L x L submatrices analogous to those in Eq. (0]l by the inverse of their top half and store the product of 
changes made by this to be multiplied with t at the end. In quantum percolation, we have to perform the 
procedure at least after each two transfer matrix multiplications (cca 10 in Anderson model), but then we 
are able to include larger samples up to L = 180. The most numerically sensitive region is the anomalous 
one, where maximum L w 100 forp > 0.1. 

A slight drawback of our method is, that even the large chosen e b still allows small conductance g 
via quantum tunneling - we can guess this value roughly by constructing a single perpendicular barrier of 
thickness 1 site, and calculate its g w 10~ n for sample sizes around L = 100 (and Ing m —25). For two 
distant barriers we get g w 10~ 21 after tunneling twice. Deep insulators must be treated carefully, in order 
to distinguish between Anderson localization and quantum tunneling. 
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A part of our finite samples creates percolating cluster of insulating sites even for ratios of these sites 
p < p c , recall that p c is the classical site percolation threshold. In classical percolation such a sample 
would have zero conductivity, in our case quantum tunneling gives small non-zero conductance of above 
mentioned orders. Neither taking these samples into account nor omitting them is fully correct. Thus 
when we calculate g in the localized regime, we prefer large samples and ratios p not quite close to p c . 
The accuracy of (In g) is somewhat spoiled and we will therefore use the less sensitive (g) and var g = 
(.9 2 ) — (s) 2 > if possible. At least this is easy to test by changing €b- If the conductance changes, we 
have tunneling or another (numerical) problem. Otherwise we have Anderson localization and the data are 
plausible. 




3 Numerical results 

At first let us describe the so called anomaly region at the center of the conduction band, i. e. for small 
energy \E\. It is known, that E — can be a special case in these models, e. g. in ID Anderson 
model the wave function gains additional symmetries, Ref. ifTTl . Fig. Q] shows the conductance as a 
function of energy, g{E), for three realizations of positions for B-type sites, with the same ratio p. We 
can see that, outside of the center, (g) slightly increases with the number of open channels (which itself 
decreases with \E\) and g oscillates with approximately constant fluctuations. The ergodic theorem states 
that these fluctuations of g are the same as those for many different samples with the same p and E, called 
statistical ensemble. We will calculate these fluctuations in the latter way. But near the band center g 
quickly decreases to some random value, much smaller than g far from center. Such behavior remains in 
localized regime, too. This anomaly is not present in the 2D Anderson model. We believe it is connected 
with the decreasing density of states (DOS) p at the band center lfTSl [181. because the ensemble averaged 
DOS (lOx scaled for convenience) p(E) is similar to (g{E)) in metallic regime, see FigQ] The lack of 
eigenstates reduces the transmission, though not in a trivial manner. We found that mean g(E;p, L) is not 
simply proportional to p(E;p, L). An analogical graph, without DOS, was published in Ref. |fl9ll , where 
even an average over 100 configurations did not remove the fluctuations completely, not to speak about the 
anomalous decay. Similar DOS, with larger disorders p, where presented in Ref. [12]. The high peak at 
E = was removed. 

Another interesting feature of the anomaly is the partial quantization of conductances, see Fig. [2] It 
was already reported in Ref. [3| for g as a function of the wave number k. Here it happens at the band 
center, for medium disorder p and it is especially well pronounced for smaller samples. The most typical 
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values of g are slightly below 2, or 4, etc. This means that we have 1, or 2 conducting channel(s) and the 
rest is almost negligible, as the peaks are highly asymmetric, practically no 2+, or 4+ conductances. For 
larger g the peaks are shifted to lower values (lower than 6, 8) and suppressed, again as in Ref. Q, where 
it was shown that the quantization disappears in the thermodynamic limit. Accordingly in our case both 
for larger \E\ and/or L this quantization is gradually smeared out. Closer inspection of samples shows that 
these conducting channels are not trivially connected to existing straight lines of conducting sites across 
the sample, because "pathways" with short detours can also give jsi2 and others can have very small g. 

There is no smooth change of P{g) with one maximum from Gaussian distribution in metal to deformed 
log-normal distribution in the insulator, as it is the case far from band center, Fig. [4] and for the Anderson 
model. Instead a part of samples enters these quantized peaks and the others enter the clearly separated 
peak of localized samples close to g — 0. With larger L more samples become localized and less samples 
remain quantized, regardless of (fixed) p. We found no critical P c (g), that would be L independent. This is 
why we are convinced that there is no MIT even at E = and the discontinuous behavior of conductances 
could have lead some authors to different statements, e. g. Ref. 0. 
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Fig. 3 Resonance peaks of g(E) for insulator. 

Further let us address the problem of resonances, manifest in the deeply localized regime in the whole 
conduction band. For smaller samples we get the well known resonance peaks with relatively high ampli- 
tudes yielding maxima close to the value g max — 2, see e. g. Ref. JTJ. 

We present g(E) semilogarithmically on a small interval of E with several typical peaks in Fig. [3] the 
maxima being roughly Lorentzians. We often get a coupled minimum - maximum pair, apparent in lng(E) 
plot. These line shapes are somewhat reminiscent of those from Azbel - Kaner [20] cyclotron resonance, e. 
g. Ref. El : just the latter are for g instead of Ing. Of course, we have no successive maxima E n = nEi, 
as we have no magnetic field and the similarity is probably occasional. There are also few isolated minima, 
even deeper than the paired ones. The minima where not reported yet, as they are not very apparent in 
linear scale. 

The values of global maxima are slightly below 2, say j« 1.95 — 1.99, in analogy to the quantization in 
the anomalous region. This means that we have one open channel, almost perfectly conducting at special 
tuned energy and the contribution of other channels can be neglected at maxima. Calculating an ensemble 
of 10 5 samples with fixed E and p up to 0.4, virtually never gives conducting samples. Nevertheless 
performing a rough energy scan of one sample and then tuning the E of chosen maximum up to 9 digits 
allows to find them already within few samples - this is a much more effective strategy for finding these 
rare events. Contrary to anomalous quantization we never see g m ax ~ 4, i. e. two channels tuned to same 
E. The peaks are homogeneously spread along the whole conduction band. 
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Now let us switch to more common strategy of statistical emnsembles. I. e. we calculate the properties 
of many samples with fixed E and p, just not too close to p c . For E = —0.53, i. e. far from band center, 
medium system sizes L and large ensembles (10 6 samples) we fitted the ratio of extended cases, Ref. [ 14|, 
r(g > 2) ex exp(— cL K ) and got the value n = 0.52(11), close to Pendry's 1/2. We performed identical 
analysis for Anderson model and got almost the same k 0.55(12). Error bars are rather large, the lnr(L) 
data get noisy for very small r and this reason disabled analogous verifications in 3D cases. Nevertheless 
these results are strong support of Pendry's conjecture k = 1/2 for both models in 2D. The g = 2 is a 
natural turning point as, apart from the anomalous band center, g > 2 usually means, that at least one or 
two channels have transmission close to 1 and g < 2 arises typically as a sum of many channels with low 
transmission; this is the well-known cause of non-analyticity of P(g), Ref. 11221 . see Fig. [4]b, c. Contrary to 
the anomalous region, the quantized cases with g w 1.99, become very rare within the ensemble approach 
for deep localization case, as the peaks in Fig. [3]are narrow. 
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Fig. 4 P(g) for E = -0.53, L = 80 and various p. 



3.1 Distributions of conductance, far from band center 

The distribution P(g) is Gaussian in metallic regime, but in crossover regime between metal and insulator, 
L w £, it differs significantly from Anderson model's distributions even outside the anomaly region. It 
strongly resembles the distributions calculated in Ref. ll23l for quasi-lD wires with rough surfaces. It is 
interesting, that assembling (still randomly) the insulating sites on one side or dispersing them in volume 
does not change typical P(g) even in the crossover regimes. We can confirm that also in our model P(g) 
can be categorized by (g) and it depends only slightly on the combination of other parameters, giving the 
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same (g), thus confirming the single parameter scaling. Just our P(g) of the same shape correspond to 
somewhat lower (g) values than theirs (Ref. [23], beware the units), but our samples are much larger. 
We recall that E = —0.53 was now chosen so that we are far from the anomalous band center and P(g) 
changes continuously with p, i. e. its maximum shifts gradually instead of being fixed to some quantized 
values as in Fig|2] 

Some noisy (because of much smaller statistics) P(g) histograms for our model were already published 
fl9l . We used ensembles of 200.000 samples. 

As usual in these models, we have non-analyticity of P(g) at g = 2. The only difference in shapes 
is, compared to wires with rough surfaces, that our P(g) in FigHH is less sharp than that in Ref. |23l : it 
turns gradually to a Gaussian and its maximum is slightly above g = 2; for p even smaller the peak shifts 
towards larger (g) . The sharper peaks in ll23l were calculated for very small number of open channels, 
which is not our case. 

3.2 Cumulants of conductance 

Let us now plot several dependences of (g) and var g on model parameters L or p. In analogy to Anderson 
model we expect that the leading term of the mean conductance g oc L in the ballistic regime (for both 
L and p very small) and g oc L° in the metallic - diffusive regime. In Fig[5]we can see that indeed for 
p w 0.1 we get a wide region where (g) is almost independent of available L, though we expect that for 
any positive p, large enough L (beyond computing possibilities) would localize all modes and there is no 
MIT. To confirm this, we calculated up to L — 160 and found a small but significant decrease, the onset of 
localization. 




In the insulating regime, there were indications, that the localization in the anomaly region could be non- 
exponential, e. g. power law or superlocalizing. It was already shown in Ref. [ 19 1 numerically, that at least 
in 2D for L large this is not the case - the (In g) vs. L plot is linear also for small \E\. We show in bottom 
part of Fig|6]that even (g(L)) in semilogarithmic scale is almost linear for larger L, thus the Anderson 
localization is present both inside (E — —0.02, symbols with broad boundary) and outside the anomaly 
region (E — —0.2, full symbols, and E = —0.53, empty ones). The deviation from linearity at smaller 
samples, see e. g. the lowest line, could be misleading. We present here also (g(L)) in the anomalous band 
center E = 0, Fig|7] It is similar to Fig|6] but the plateau-like regions, i. e. (quasi)metallic parts around 
maxima are shorter and shifted to lower disorders p and higher (g). We can see exponential decay, i. e. 
Anderson localization, for p > 1 and thus no MIT for p = 0.26 and E = 0, reported in Ref. 0. 
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We will further analyze the conductance fluctuations. In metallic regime the half width of P(g) Gaus- 
sians is an universal constant, depending only on physical symmetries, boundary conditions (BC) and 
dimensionality, but independent of energy, system size and disorder - for medium p and L; in higher di- 
mensionalities d > 2 even for the largest L available. We use the hard wall BC in perpendicular direction, 
quantizing the last term in Eq. (O. The slow dependence of y/var g on L was already published [ 19], we 
prefer var g(p), where the theoretical prediction approaches the index point, see FigJH] 




Fig. 7 Mean conductance as a function of system 
size in the band center. 



CO 

> 



1.2 
1.0 

'0.8 
0.6 
0.4 
0.2 



ft E=-0.53 






o-o L=64 : 




EH3 L=80 




<^>L=100 Z 




^L=120 - 







0.05 



0.1 



0.15 



0.2 



0.25 



Fig. 8 Conductance fluctuations vs. disorder for 
various L. Broken line is theoretical UCF. 



The universal conductance fluctuations (UCF) predicted for 2D hard wall case ought to be var g w 
4 * 0.1856 = 0.7424 in our units, e. g. Ref. EHEU. One can see in FigH]that after the ballistic peak 
there is a region, where var g only very slightly decreases with p and this happens around the theoretically 
predicted value. This graph is very similar to the one calculated for Anderson model 11261 . This is true only 
outside the anomaly region; var g at maxima in Fig. [7] become non-universal and at least by one order 
of magnitude greater than UCF - in band center we are still in ballistic regime. Nevertheless, the metallic 
regime is not just the only one exactly solved yet, predicting Gaussian P(g) with known universal halfwidth 
and (g) oc L d ~ 2 ; it is also the most robust regime for the mutually related models with disorder-induced 
localization. 



4 Conclusions 

We have shown, that the transfer matrix method suits well for studies of the 2D quantum percolation model 
as a limit of binary alloy. The regions of interest can be divided according to following two criteria. The 
first one is the regime of transport, given mainly by disorder p: ballistic, quasi-metallic, crossover, localized 
and deeply localized regimes. Another one is according to energies E: close to band center (small \E\, 
called anomalous region) and far from it. 

The statements of MIT present in 2D quantum percolation model, Refs. Q and [6 ], were achieved with 
ID linear leads and smaller conduction band than ours. We have leads of thickness L > 1. It was shown in 
Ref. [27] that narrow leads give different results, at least in 2D Anderson model, than the wide ones. On 
contrary in Ref. |4] the absence of MIT was shown for E — 0.5. We did not find the MIT either, even at 
E = 0. Increasing g(L), that would give rise to positive values of the scaling function (3(g) = d\ng/d\n.L 
(and MIT at (3(g c ) = 0), appears to be just finite-size effect. We believe that (3(g) is always negative|4|. 
The anomaly region near the center of the conduction band retains Anderson localization in large samples 
and we found neither sublocalization (power-law decay of g(L)), nor superlocalization (over-exponential) 
for large L at any E, p. Instead, deviations appear that are again just finite-size effects, unusually magnified 
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by the novel partial quantization of g. It is a question, whether we should even call the conducting regime 
at band center (quasi)metallic, as the plateau of mean g(L) is too small, no UCF are observed and the onset 
of exponential decay (localization) appears already for g > 2. The only reason left is that the peaks of 
P(g) apart from quantization region, say those with (g) > 10, are almost Gaussian. Anomalous behavior, 
seen theoretically in ID at E = and sometimes conjectured to higher dimensionalities, is not found 
numerically for Anderson model in 2D and 3D, but it appears in the quantum percolation model for small 
and medium L. We suppose that the decay of g(E) at band center is also connected with the decrease of 
DOS in the sense that effectively non-zero p is a necessary condition for non-zero g. 

Most Lorentzian (or Cauchy) resonance peaks in the strongly localized case are of comparable relative 
height in semilog-scale and thus they do not make the samples always conductive in our calculations of 
deep insulators, i. e. for both p and L large. Pairs of minimum and maximum - or vice versa - occur in this 
scale that were not reported yet. 

The probability r(g > 2) oc exp(— cL K ) of rare delocalized samples confirms Pendry's conjecture |[T4l 
k = 1/2 both for 2D quantum percolation far from band center and for 2D Anderson model. This is the 
first numerical confirmation of the conjecture. The rare, i. e. conducting cases within deep insulators, 
close to classical percolation threshold, can be found by tuning the energy at maxima in g(E) plots even 
for medium sized samples. 

The distributions of conductance in the crossover regime far from band center are very similar to those of 
quasi- ID wires with rough surfaces and they are different from those of Anderson model. In the metallic 
phase the mean conductance is practically L independent in a large region and, except band center, its 
variance is close to the universal (UCF) value predicted theoretically 11241 . 

Concludingly let us compare the Anderson model and the quantum percolation one. We are convinced 
that they belong to the same universality class, i. e. they undergo no MIT for d < 2. Running forward, 
in our most recent paper l28l we show that the models have identical critical exponent in 3D MIT. Back 
in 2D, far from band center they behave very similarly, apart of tiny differences in crossover regime. At 
band center, quantum percolation model has very strong finite-size effects, somewhat questioning even the 
single parameter scaling hypothesis. Such effects are not present in Anderson model. 

This work was supported by Grant VEGA Nr. 2/6069/26. 
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